Understanding physical climate risks and their implication for community adaptation in the borana zone of southern Ethiopia using mixed-methods research

Pastoralists in the Borana zone of southern Ethiopia are grappling with climatic risks and impacts, and their livelihoods appear to be precarious. The long-term changes in precipitation, temperature and extreme events are under-researched aspects of climate risk in the study area. Climate observations, local people’s experiences, and memories of weather and drought patterns were analysed in this work to better understand the climate risks in the study region and provide actionable knowledge to facilitate adaptations. In southern Ethiopia, the primary drivers of climate risks are rainfall variabilities, rising temperatures, and drought. The annual rainfall variability observed in the study area (20–35%) is greater than the overall estimate for Eastern Africa (15–25%), indicative of a heightened climate risk in this area. Furthermore, seasonal and intra-seasonal rainfall variabilities are being aggravated by rising temperatures, leading to increased frequency and magnitude of extreme events, particularly in the lowland area. These changes, specifically droughts of different intensities occurring every 2–3 years, erode the livelihood of the pastoralists. The lack of consensus among stakeholders as to the causes and aggravating factors of these climate changes impedes adaptation actions. The authors recommend the importance of initiating a participatory platform that will facilitate discussion on climate change to create a common understanding of the problem and relate it to public policy. The use of gridded precipitation and temperature data is recommended in future research to improve the availability of climate information for climate risk management. In addition, the use of mixed methods and local knowledge, as demonstrated in this study, is highly recommended to better understand climate risks, particularly when there is a limited availability and quality of long-time-series climate data.

Stations that record precipitation data are classified as class 4, while those that record both precipitation and temperature data are considered class 3. Class 1-2 s station records relative humidity and wind direction/speed data in addition to precipitation and temperature data. Most of the stations in the study area are considered class 4 or 3. There is only one Class 1-2 station included. All available station-recorded daily data used in this study, spanning from 1986 to 2018, were collected from the NMA. Figure 1 summarizes the the data availability for each station.
Remotely sensed data. Moderate-Resolution Imaging Spectro-radiometer (eMODIS) normalized difference vegetation index (NDVI) data processed by the United States Geological Survey (USGS) were widely used by the Famine Early Warning Systems Network (FEWSNET) for drought monitoring. The Dekadal (10 days) raster datasets from 2002 to 2020 were acquired from the FEWSNET data portal to ascertain trends in NDVI and complement the observed data.
Household data. Household demographic characteristics, common livelihood risks in order of importance, and perceptions of drought, climate variability, and change trend, among other questions, were included in the closed-ended survey instrument. Rainfall (amount, rainy days, timing-onset or cessation, and distribution), temperature, drought (frequency, duration, severity, and intensity), and flood (frequency) were key parameters included in the instrument. The past 10 years were used as the "recall period" for the trends of the perception questions given the age variation of the households. The instrument was digitized into Open Data Kit (ODK) tools and facilitated mobile data collection using a smartphone. The detailed sampling procedure, sample size, and sample distribution in the study area are presented in Annex 1.
Qualitative data. Focus group discussions (FGDs) and key informant interviews (KIIs) were the qualitative assessment methods applied in this study to complement the climate and household data. A mixed focus group interview was facilitated due to the extensive nature of the study area. Women, men, youth, elders, and various ethnic group representatives constituted the FGD. Community sketch/resources map, common livelihood risks ranked in order of importance, and perceptions of drought, climate variability, and change trends, among other factors, were included among the semi structured questions. The past 30 years were used as the reference period for the perception questions. Twelve FGDs were facilitated in total, one per district, with the exception of Teletel (two FGDs) because of its size and diversity of livelihoods (refer to Annex 2 for a map of the FGD distribution). Diverse groups endowed with first-hand knowledge of the subject matter were used as key informants. Traditional grazing managers locally called "Jarsa Deheda" were one of the categories of KIs. The research team interviewed four of the five recognized grazing managers in Borana, i.e., Globo, Malebee, Wayama, and Gomole.  www.nature.com/scientificreports/ The Dire grazing manager was not available due to his medical condition. Some of the indigenous weather forecast practitioners were also used as key informants. One of the known animal behaviourists ("Waragu") and cyclical calendar practitioners in Borana ("Mara") were interviewed as KIs. Discussions with representatives of relevant government departments of the Borana zone were further used as KIIs. A total of seven KIIs were facilitated by the research team. Semistructured questions similar to the FGDs were included to ascertain the key informants' understanding. Drought (frequency, duration, severity, and intensity), flood (frequency), rainfall (amount, rainy days, timingonset or cessation, and distribution) and temperature were key parameters included in the instrument. Likewise to the focus group discussions, the past 30 years were used as the reference period for the perception questions. The "timeline" of major drought events was expanded to 60 years by most of the KIs during the interview process.
Secondary data sources. The United Nations Office for Disaster Risk Reduction (UNDRR) has collated disasterrelated losses and damages that prevailed in Ethiopia from 1994 to 2013 from various sources. These digitized data are available on the Humanitarian Data Exchange (HDX) platform. Major historical disaster data relevant to the study were downloaded from HDX. These data were used to complement the historical timeline of the key informants to characterize prevailing disasters in the study.
Quality control of meteorological data. The acquired climate data were further assessed for completeness and quality before analysis. All of the stations suffered from data gaps and incomplete historical records, with some lapses in reporting being larger than others. Figure 2 shows the total number of years that each station has reported data relative to the number years of complete observations (no missing data). The annual observations of a given station were considered "complete" when rainfall or temperature data were available for 10 months or more in a given year.
Climate data often have missing information, which limits their use 31 . Stations with significant years of incomplete data are believed to affect the overall statistical results obtained using the time-series data. Estimating missing data is a routine procedure to address these gaps, and there are standard methods of doing so. The arithmetic mean method, normal ratio method, modified normal ratio method, and interpolation method are some of the most widely practised methods. Missing rainfall can be determined using a simple arithmetic average if the normal annual precipitation at various stations is within 10% of the normal precipitation 32 . If the normal precipitation varies considerably (exceeding 10%), then missing rainfall can be estimated by weighting the precipitation at various stations by the ratios of the normal annual precipitation 33 . The other method used to fill missing meteorological variables is the use of datasets from other selected stations in the surrounding area and the application of appropriate spatial interpolation methods 34 . Most of these methods, however, require the availability of well-established climatological normal information of base stations as well as neighbouring stations with grossly missing data in the study area.
In addition, a considerable estimation level is required to fill in the prevailing station data gaps that affect parameter estimations. This study used different approaches to address the data gaps in the absence of effective estimations. First, we prioritized stations with a reasonably complete data. For precipitation, these stations were Chewe Bet, Finchwoha, Yabello, and Moyale. For temperature, only one station has acceptable data coverage for  www.nature.com/scientificreports/ prioritization: Moyale. Second, we decided on a common window period for the analysis of prioritized stations. The 1990-2017 period had a lesser degree of missing values and thus a smaller estimation requirement than other periods across the stations. Hence, this period was chosen as the common window for the analysis of prioritized stations. the missing monthly rainfall data were filled using the corresponding monthly median of each station. The median was selected as it is a better representation of the central tendency, relative to the arithmetic mean, given the prevailing extreme values of rainfall data in the study area. The double mass curve (DMC) is the most commonly used graphical approach used to test the homogeneity of time-series data at prioritized stations, as presented in Fig. 3. Its low data requirement makes this technique popular, and the method can indicate a change in the slope of data as a breakpoint from which the data exhibit inconsistency. The prioritized stations were found to be consistent according to the double mass curve analysis. Some of the minor breakpoints observed at some of the stations along the trajectory emanated from the data estimations or from changes in the observing conditions, practices, or procedures over time.
Data analysis methods. Climate data analysis. Descriptive statistics are a powerful tool for summarizing different climatic parameters within specific time periods. The mean and median are a useful measures of the central value of observed elements such as precipitation and temperature and is used to approximate normality. The standard deviation (SD), standard error (SE), and coefficient of variation (CV) are among the selected dispersion statistics used to depict the extent of climate variability compared to the normal state. The selected statistics were calculated for a series of different observational periods (i.e., the monthly, seasonal and annual periods).
Parametric trend analysis (statistical or graphical) methods are used to ascertain trends in observed and remotely sensed time-series data. The nonparametric Mann-Kendall (M-K) test and Sen's slope estimator were analysed herein using XLSTAT 2017 software to determine the direction, statistical significance, and magnitude of the trends of the observed and remotely sensed data.
The rainfall anomaly index (RAI) was used to characterize the severity of drought stress in the study area. The RAI is a simple drought index used to analyse the frequency and intensity of dry and wet seasons or years and was developed by van Rooy in 1965 32 . It incorporates a ranking procedure to assign magnitudes to positive and negative precipitation anomalies. The equation of the index is shown as follows: Seasonal/annual rainfall anomaly index (RAI): where: N : current seasonal/yearly rainfall for the RAI generation period, N : seasonal/yearly average rainfall in the historical series (mm), M : average of the ten highest seasonal/yearly precipitation values in the historical series (mm), X : average of the ten lowest seasonal/yearly precipitation values in the historical series (mm). RAI analysis can further help researchers determine the probability of occurrence of dry conditions, including the drought and recurrence intervals, in the study area. The probability of drought occurrence by severity was calculated using the equation below.
Probability of occurrence of dry conditions: where: P: absolute probability of occurrence (severity), n: number of event occurrences (RAI < − 0.99) of drought, N: Total number of possible occurrences (the study period).
(1) www.nature.com/scientificreports/ The dry and drought condition recurrence intervals or return periods for each of the stations were computed from the RAI based on the formula below: Recurrence interval or return period: Household and qualitative data analysis. The frequency distribution and simple graphic analysis were used to analyse the household data using a statistical package for the social sciences (SPSS) and Microsoft Excel. Each of the focus group discussions and key informant interview audio records was transcribed as the first step of the qualitative analysis. The transcriptions were further edited, synthesized, and analysed using the content analysis technique. The mapping, ranking, timelines, trend analysis, and ratings collated from the participatory tools provided perspectives for the subsequent qualitative data analysis.

Results
Rainfall and temperature descriptive analyses. The mean, median and SE statistics for annual and seasonal precipitation were summarized for the selected stations and are presented in Table 1. The contribution of each rainy season to the annual rainfall total was also calculated and is provided in Table. In addition, the annual and seasonal average temperatures for the Moyale station are presented in Table 1.
The coefficient of variation was calculated and presented in Table 2 to show the year-to-year or season-toseason rainfall variability. These statistics consider the deviations from the averages regardless of whether a station experienced high or low rainfall and can be used for relative comparisons of variabilities among stations. According to a cited work 35 , the CV can be used to classify the degree of variability of rainfall events as less (CV < 20), moderate (20 < CV < 30), high (CV > 30), very high (CV > 40%) and extremely high (CV > 70%) interannual rainfall variability.
Rainfall and temperature trend analyses. Parametric trend analysis. The total annual and seasonal rainfall totals of a given station are plotted in a line graph over time to ascertain any observable trends. A similar analysis was performed for the Moyale station mean temperatures. A simple linear regression equation was generated using MS Excel to predict the annual precipitation in mm or mean temperature in • C for a unit change of the dependent variable (time in a year). In addition, linear trend lines and R-squared statistics were generated Time(Years) Number of events(RAI less than − 0.49 occurred) Table 1. Summary statistics of temperature, precipitation and contribution of each season to annual rainfall (1990-2017). Source (Author): Computed based on data from the NMA of Ethiopia. Long rainy season, "Genna" (Mar-May); short rainy season, "Hagaya" (Sep-Nov); cold dry season, "Adolessa" (June-Aug) and long dry season, "Bonna Hagaya" (Dec-Feb). a Median. www.nature.com/scientificreports/ from the built-in tool to determine the direction of change and accuracy of the fit. The slope of the equation (trend) and R-square value for each station were extracted from line graphs and are presented in Table 3.
Nonparametric trend analysis. Trend analyses were conducted using MK and Sen's slope estimators for seasonal and annual precipitation at each station. The summary statistics obtained for the main, short-season, and annual precipitation are presented in Table 4. Similarly, the MK trend analysis results for the maximum and minimum temperatures at the Moyale station, including the monthly trends, are presented in Tables 5, 6.
In the Mann-Kendall test results (S), a positive S value indicates an increasing trend, while a negative value indicates a downwards trend based on its significance. The significance test is based on rejecting or not rejecting the null hypothesis (Ho, i.e., no trend in the data series) at a 5% significance level (alpha value). If the p value is below the threshold of significance (p < 0.05), then the null hypothesis can be rejected, and the prevailing trend (increasing or decreasing) can be considered significant. Kendall's tau is a nonparametric rank correlation coefficient that measures the strength of the relationship between the observed element (e.g., precipitation or temperature) and time. The higher the value of the coefficient is, the stronger the relationship is, and vice versa. The estimated Sen's slope further provides further information about the trend of the observed element. A positive Sen's slope reveals an upwards trend, while a negative Sen's slope suggests a downwards trend.
Characterizing drought and its severity. The seasonal rainfall anomaly index (ARAI) was calculated based on Eq. (1) presented above for each of the two main rainy seasons of the observed year at selected stations. The index values were further assessed against the 9-class scheme proposed by Van Rooy 36 , ranging from extremely wet (≥ 3.00), very wet (2.00 to 2.99), moderately wet (1.00 to 1.99), and slightly wet (0.50 to 0.99) to near normal (− 0.49 to 0.49), slightly dry (− 0.99 to − 0.50), moderately dry (− 1.99 to − 1.00), very dry (− 2.99 to − 2.00) and extremely dry (≤ − 3.00). Counts of wet and dry rainfall conditions corresponding to various severities for the two seasons are summarized in Table 7.
The overall observed years at each station were classified into three categories, dry, normal, and wet years, based on the ARAI. In addition, the probability of occurrence of each category was calculated based on Eq. (2). The drought recurrence interval or return period was computed using Eq. (3). The result is presented in Table 8.  Fig. 4. A simple linear regression equation was generated to predict the monthly NDVI for a unit change of the dependent variable (time in a year). In addition, linear trend lines and R-squared statistics were generated from the built-in tool to determine the direction of change and accuracy of the fit. This information corresponding to the annual period and two seasons is incorporated in the line graphs. The nonparametric trend analysis method incorporating the M-K test and Sen's slope was used here to determine the direction, statistical significance and magnitude of the trends of the Dekadal NDVI, and the results are presented in Table 9.     Fig. 5. The youngest to oldest age groups (20 to 95 years old) were interviewed to help bridge views and opinions across generational gaps. Most of the respondents (78%) had resided in the area for an extended period (20 years was a modal value) in the same location, which helped us reasonably assume their familiarity with observable changes in their area. Only 12% of the respondents claimed to have migrated from other areas in the last 10 years or less. Households clearly indicated that the frequency and severity of drought and temperature have increased compared to the last 10 years. Households were also undivided in their responses regarding the reduction of the rainfall amount and number of days in the two main seasons. They also unanimously agreed on the erratic nature of the rainfall pattern in recent time compared to over the past 10 years. The frequency of flood events and extreme temperatures were overall perceived to have increased, although the response of the households were not as united as those corresponding to the previous indicators.
The KII assessment provided further insight into climate change parameters in addition to affirming the household analysis results presented in Fig. 5. The informants widely cited rainfall pattern changes as evidence of their observed change. The informants compared rainfall patterns in terms of their timing, amounts, and distributions in the present day and over the past 30 years to substantiate their claims, as presented in Table 10.
The informants were unanimous in terms of the rising temperature trend in the area and linked this trend with the changing pattern of rainfall. In the past, the weather during the long dry season was often characterized by relatively high temperatures during the day and low temperatures, including an uncomfortable coldness, at night. However, low temperatures (coldness) are observed during both day and night in the cold dry season. In recent times, however, increasing temperatures have been observed regardless of the season or day/night.
Drought recurrence and severity were other parameters the informants compared to support their assertion of climate change. For example, Malbee grazing managers asserted that as droughts became more frequent, even before livestock keepers recovered from one drought, another drought hit the area. The historical timeline collated from the key informant interviews provided further insights into the temporal dimensions of major hazards. The synthesis of major drought events and associated disasters the population has experienced since the early 1960s and their chronological characteristics corresponding to the Gadaa calendar are presented in Table 9. M-K test of NDVI anomalies from 2002 to 2020 in the Borena zone. www.nature.com/scientificreports/ Annex 5. This synthesis can help us understand the causes and effects of the events and put the research into a historical perspective. Drought and dry spells are among the major hazards threatening the livelihoods of the population group. The sample household ranking order of these common hazards further augments key informant assessments. Erratic rainfall amounts and distributions, including drought events, ranked first among the common hazards faced by sample households. Our analysis of extracted digitized data recorded from 1994 to 2013 revealed that drought was the most common disaster during this period. Drought events were widely distributed spatially across the studied zone, while other hazards affected localized areas of the zone, as presented in Annex 6.
Informants were quite divided as to the causes and contributing factors of climate change. The majority of the elders, however, strongly believed that climate change is an "act of God" and did not subscribe to the scientific basis of human-interference (anthropogenic) causes of climate change. Some of the informants attributed the changes to the weakening of traditional customs and practices.
One of the grazing managers and most of the KIs from the government acknowledged the link between deforestation and climate change. The KIs highlighted the gradual deforestation pocket of the forested area in Borana (Miyo, Arero, and Yabello) due to various human activities and believed that these changes contributed to local climate change, particularly to the warming of temperatures. Land use and land cover change (LULC) observed in the area due to internal and external factors might also have contributed to or aggravated local climate change.
The informants from government institutions did not discount the aspect of the global climate change process, i.e., greenhouse gas (GHG) emissions aggravating the local climate conditions. Methane gas (CH4) emitted through livestock disposal has its own share of concern given that livestock-keeping is the predominant means of maintaining livelihood in the area.

Discussion
High rainfall variabilities observed in the study area are similar to the patterns recorded in Eastern African drylands, which have very high inter-annual variability in the length of the growing season, particularly for areas with bimodal distributions 37 . East Africa experiences some of the largest inter-annual rainfall variations in the world, and there are many different drivers of this variability 38 . However, the higher temporal and spatial rainfall coefficients of variation in the study area compared to the estimate (15-25%) for Eastern Africa 37 are concerning and indicative of the increasing extreme rainfall events.
No statistically significant trends were observed in the seasonal, or annual mean rainfall at the four analysed stations. However, the observed downwards trend of the main-rainy-season precipitation among the four stations mirrored the widely held perception of the community. The upwards trend of the short rain season observed, although not statistically significant, was quite contrary to the household responses and qualitative and NDVI trend analyses. For example, the vegetation condition trend (NDVI) of the short rainy season showed a statistically significant downwards trend that was strongly related to the rainfall received in semiarid regions of Africa 39 . Qualitative and household analyses further conclusively suggested that the amount of rainfall decreased and became insufficient, notably for the short rainy season, compared to the past. Some studies in Ethiopia have identified downwards trends in precipitation in some parts of the country, mainly in the eastern, southern, and www.nature.com/scientificreports/ south-eastern regions 40,41 . The key informant comparisons of past and current rainfall pattern timing, amounts, and distributions presented in Table 10 provide substantial and contextual information as to unfavourable precipitation, temperature, and drought trends. Quantifying observed trends in precipitation is more challenging than quantifying those in temperature due to the lack of reliable data and variability on multiple temporal and spatial scales. However, significant declines in the Long Rains from March to May over the Horn of Africa were observed in the 1980s and 1990s 42 , followed by a recent recovery 43 . The trend analyses of the mean monthly, seasonal, and annual maximum and minimum temperatures at Moyale station showed increasing (warming) trends overall, although the results were mixed in terms of their statistical significance. The February, March, May, June, and July monthly maximum mean temperatures were statistically significant. Similarly, the minimum mean temperature for the June-October period exhibited positive and statistically significant trends. Seasonally, statistically significant warming of maximum temperatures was observed in the cold and long dry seasons, while a similar pattern was observed in the cold dry season and short rainy season for the minimum temperature at the 5% level of significance. The key informants were highly unanimous about the warming of temperatures in the study area regardless of the season or time of day/night in recent times. Households also clearly indicated that temperatures, including extreme temperatures, have increased in the area compared to the past 10 years. Early analyses 44 suggested that the temperature increases in the semiarid lowland regions of Ethiopia were more pronounced than the national average changes as well as those in highland regions in the country. These regions, particularly Borana, Guji, and South Omo, experienced temperature increases of 0.40 °C per decade from 1950 to 2000, further supporting the findings of this study.
The pattern observed in this study is in conformity with the two broader consensuses and evidence of climate change in Eastern Africa 45 . Increasing temperature trends from climate baselines are among the broad pieces of evidence obtained at the regional level. Second, the overall rainfall trend assessments were weak and hard to detect despite either increasing or increasing trends observed depending on the specific locations and seasons.
The observed high rainfall variability entails a high degree of rainy-season uncertainty in the study area in general, particularly for the short rainy season. This finding implies that the study area, particularly the low-land area, is exposed to a relatively high degree of climate risk, particularly due to extreme precipitation, because of the relatively high inter-annual and seasonal rainfall variabilities and relatively low precipitation. Our seasonal rainfall anomaly analysis suggests that the study area faced probabilities of experiencing dry conditions at different intensities (− 0.5 and less) during the main and short rainy seasons of 35-50% and 35-65%, respectively. This probability implies that 4-5 and 4-6 dry events of different intensities occurred every 10 years during the main and short rainy seasons, respectively. The stations located in lowlands had higher dry-year frequencies (including higher intensities) than the stations in mid/higher elevations, linked with their exposure to a relatively low amount of seasonal rainfall and high variabilities. On average, droughts of different intensity degrees occur every 2 to 3 years in consecutive or non-consecutive seasons in the study area, as the recurrence interval or return period calculation results of the observed data suggest. The observed temperature increases will aggravate the prevailing dry conditions by enhancing evaporation, thereby reducing surface water and drying out the soils and vegetation. Rainfall and evapotranspiration are two major climatic factors affecting agricultural production 46 , and agricultural water resources face two major problems. Climate extremes significantly impact livestock productivity in eastern and western Africa. The increasing frequency and intensity of droughts, changes in water availability, and increasing patterns of temperature and rainfall variability all profoundly threaten the livelihoods of people living in drought-prone areas and the existence of society in arid and semiarid remote regions 47,48 .
The empirical results reported herein should be considered in light of some limitations. The climate risk analysis presented in this paper was based on limited observed rainfall and temperature data due to the unavailability of high-quality data from meteorological stations in the area of interest. In this research, we performed the necessary quality and homogeneity tests to ensure the integrity of the results. The use of gridded data that blend station-observed data and satellite observations or reanalysis data is recommended in a future study to broaden the scope and application of the findings of this research.

Conclusion and recommendation
Recently, unfavourable drought, climate variability, and climate change trends have been observed in the study area, according to most households and key informants. Climate data provided further credence to these household and community perceptions regarding warming temperature trends, increasing rainfall variability, and recurring dry and drought conditions. However, the climate data results did not show any unified pattern regarding the local characterization due to the limited availability and quality of long-time-series climate data in the study area. These findings highlight the need to integrate local and indigenous knowledge in climate risk studies in addition to measured climate data.
There was no consensus among the stakeholders on the causes and aggravating factors of the observed changes, though these factors have implications for adaptation actions. The authors recommend the importance of participatory workshops on climate change and its impact on the Borana pastoral system by bringing actors and the affected populations together and forging paths forward. Such workshops will create a common understanding of the problem (framing the problem) related to public policy. The workshops will further provide opportunities to create platforms hitherto absent in the area that can regularly meet to discuss the issue of climate change. The analysis and synthesis of the observed climate and qualitative data highlighted in this study serve as valuable input for building consensus.
Drought and dry spells were the major hazards cyclically grappled with by the pastoralists and agro-pastoralists in the study area. Drought events are widely distributed spatially across the studied zone compared to other hazards. A historical timeline collated from the key informants and focus group discussions provided further insights into the temporal dimensions of major hazards. The analysis captured major drought events www.nature.com/scientificreports/ the population has experienced since the early 1960s and their characteristics. The major drought events clearly demonstrate the "downwards spiral" of the livelihoods of the population. The RAI analysis suggested that drought events with different intensity degrees are recurring in 2-to 4-year intervals in the study area 18 , reiterating concern for the vulnerability of pastoral populations in the arid and semiarid zones of East Africa to recurring drought events despite the wide range of socioeconomic and cultural benefits from livestock rearing. This research used drought indices that quantify climatic anomalies from observational and remotely sensed data to characterize the drought severity, duration and frequency; to date, such an analysis has not been performed. One cited study 49 found indices to be a very useful tool in communicating complex climatic functions, including drought severity, in an easily understandable way to a wider audience. This study demonstrates the added value of such tools in decision-making in the study area. This kind of information communication tool is urgently required as the frequency, intensity, and cost of natural hazards are expected to further increase in the coming decade, as indicated by the IPCC in its fifth assessment report 50 .